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Abstract . The aeroassisted flight experiment (AFE) refers to 
a spacecraft to be launched and then recovered by the space 
shuttle in 1994. It simulates a transfer from a geosynchronous 
Earth orbit (GEO) to a low Earth orbit (LEO). Specifically, the 
AFE spacecraft is released from the space shuttle and is 
accelerated by means of a solid rocket motor toward Earth, so as 
to achieve atmospheric entry conditions close to those of a 
spacecraft returning from GEO. Following the atmospheric pass, 
the AFE spacecraft ascends to the specified LEO via an 
intermediate parking Earth orbit (PEO) . The final maneuver 
includes the rendezvous with and the capture by the space 
shuttle. The entry and exit orbital planes of the AFE spacecraft 
are identical with the orbital plane of the space shuttle. 

In this report, with reference to the AFE spacecraft, an 
actual GEO-to-LEO transfer is considered and optimal trajectories 
are determined by minimizing the total characteristic velocity. 
The optimization is performed with respect to the time history of 
the controls (angle of attack and angle of bank) , the entry path 
inclination and the flight time being free. Two transfer 
maneuvers are considered: (DA) direct ascent to LEO; (IA) 

indirect ascent to LEO via PEO. 

While the motion of the AFE spacecraft in a 3D-space is 
described by a system of six ODEs, substantial simplifications 
aj -0 possible if one exploits these facts: (i) the instantaneous 

orbital plane is nearly identical with the initial orbital plane; 
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(ii) the bank angle is small; and (iii) the Earth's angular 
velocity is relatively small. Under these assumptions, the 
complete system can be decoupled into two subsystems, one 
describing the longitudinal motion and one describing the lateral 
motion. 

The angle of attack history, the entry path inclination, and 
the flight time are determined via the longitudinal motion 
subsystem; in this subsystem, the total characteristic velocity 
is minimized subject to the specified LEO requirement. The angle 
of bank history is determined via the lateral motion subsystem; 
in this subsystem, the difference between the instantaneous bank 
angle and a constant bank angle is minimized in the least square 
sense subject to the specified orbital inclination requirement. 

It is shown that both the angle of attack and the angle of 
bank are constant. This result has considerable importance in the 
design of nominal trajectories to be used in the guidance of AFE 
and AOT vehicles. 

Key Words . Flight mechanics, astrodynamics , hypervelocity 
flight, aeroassisted orbital transfer, aeroassisted flight experiment, 
optimal trajectories, guidance trajectories, decomposition techniques, 
longitudinal motion, lateral motion, sequential gradient- 
restoration algorithm, nonlinear two-point boundary-value problems. 
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Notations 

C D = drag coefficient; 

C L = lift coefficient; 

C LP = projected lift coefficient; 

= moment coefficient; 

D = drag, N; 

DP = dynamic pressure, N/m 2 ; 

E = lift-to-drag ratio modulus; 

g = local acceleration of gravity, m/sec 2 ; 

h = altitude, m; 

h a = thickness of the atmosphere, m; 

HR = heating rate, W/m 2 ; 
i = orbital inclination, rad; 

L = lift, N; 
m = mass, kg; 

r = radial distance from the center of the Earth, m; 
r e = radius of the Earth, m; 

r = = radius of the outer edge of the atmosphere, m; 

CL 

o 

S = reference surface area, m ,• 
t = T/t = dimensionless time; 

T = running time, sec; 

V = velocity, m/sec; 

V a = circular velocity at r = r a , m/sec; 

V* = reference velocity, m/s ec; 

a = angle of attack, rad; 

y = path inclination, rad; 

r| = wedge angle, rad; 
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9 = longitude, rad; 

y = bank angle, rad; 

y e = Earth's gravitational constant, m 3 /sec 2 ; 

p = air density, kg/m 3 ; 

o = reference air density, kg/m 3 ; 

x = final time, sec; 

$ = latitude, rad; 

X = heading angle, rad; 

go = angular velocity of the Earth, rad/sec; 

q = longitude of the ascending node, rad; 

AC l = lift coefficient range; 

AC Lp = projected lift coefficient range; 

AV = characteristic velocity, m/sec. 

Subscripts 

0 = entry into the atmosphere; 

1 = exit from the atmosphere; 

00 = exit from the initial orbit; 

11 = entry into the final orbit; 

22 = entry into the parking orbit. 

Superscripts 

• = derivative with respect to dimensionless time; 

~ = variable computed in an inertial system. 

Acronyms 

AFE = aeroassisted flight experiment; 

AOT = aeroassisted orbital transfer; 

DA = direct ascent; 
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DAOT = direct ascent optimal trajectory; 

GEO = geosynchronous Earth orbit; 

HEO = high Earth orbit; 

IA = indirect ascent; 

IAOT = indirect ascent optimal trajectory; 

IART = indirect ascent reference trajectory; 

LEO = low Earth orbit; 

ODE = ordinary differential equation; 

OT = optimal trajectory; 

PEO = parking Earth orbit; 

RT = reference trajectory; 

SGRA = sequential gradient-restoration algorithm; 
TPBVP = two-point boundary-value problem. 
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1 . Introduction 

The field of aeroassisted orbital transfer (AOT) has 
received considerable attention in recent years. For a spacecraft 
which must be transferred from a high Earth orbit (HEO) to a low 
Earth orbit (LEO) , considerable savings in characteristic 
velocity are possible if the AOT mode is employed instead of the 
all-propulsive mode. In turn, this leads to propellant 
savings/payload increases. 

In the all -propulsive mode or Hohmann transfer mode, the 
spacecraft navigates in the region between HEO and LEO. It must 
be braked propulsively twice, once to deorbit from HEO and once 
to achieve circularization into LEO. 

In the AOT mode or synergistic mode, the spacecraft deorbits 
from HEO, undershoots LEO, enters the Earth's atmosphere, and 
uses the aerodynamic forces in order to deplete excess velocity. 
Then, after exiting the atmosphere, it reaches LEO with 
nearly-circular velocity. The resulting savings in propellant 
weight are beneficial; however, added thermal protection is 
needed to cope with the heating rates generated while the 
spacecraft traverses the upper atmosphere at hypervelocity 
speeds . 

Aeroassisted orbital transfer is not only important for 
HEO-to-LEO transfer maneuvers, but may prove to be indispensable 
to future planetary flights. In particular, this statement refers 
to lunar return vehicles, Mars exploration vehicles, and Mars 
return vehicles. Indeed it is known that, for a round-trip 
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Earth-to-Mars mission, the total characteristic velocity of the 
AOT mode is about half that of the all-propulsive mode. 

While the AOT prospects are clearly bright (Ref. 1), to take 
proper advantage of them it is necessary that guidance and 
control systems be designed with care. In turn, this requires the 
previous study of optimal trajectories, since they supply the 
ideal benchmark that guidance trajectories should strive to 
approach (Refs. 2-7). 

1.1. Aeroassisted Flight Experiment . Because the AOT idea is 
yet untested, NASA has planned an aeroassisted flight experiment 
(AFE) involving a spacecraft to be launched and then recovered by 
the space shuttle (Refs. 8-9). The experiment simulates a 
transfer from a geosynchronous Earth orbit (GEO) to a low Earth 
orbit (LEO) and is tentatively scheduled for 1994. Specifically, 
the AFE spacecraft is released from the space shuttle and is 
accelerated by means of a solid rocket motor toward Earth, so as 
to achieve atmospheric entry conditions close to those of a 
spacecraft returning from GEO. During the atmospheric pass, 
aerodynamic/ thermal data are gathered for use in designing AOT 
vehicles. Following the atmospheric pass, the AFE spacecraft 
ascends to the specified LEO via an intermediate parking Earth 
orbit (PEO) . The final maneuver includes the rendezvous with and 
the capture by the space shuttle. Clearly, the entry and exit 
orbital planes of the AFE spacecraft are identical with the 
orbital plane of the space shuttle. 
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The configuration chosen for the AFE spacecraft is the so 
called raked-cone configuration (Fig. 1) , which is dominated by 
heating rate considerations racher than by aerodynamic performance 
considerations. Therefore, the lift-to-drag ratio modulus is very 
low, E = 0.28, at the design angle of attack, a = 17 deg; 
also, the lift coefficient range is quite low, -0.47 < C L < -0.21, 
in the angle of attack range 2 7 >_ a >_ 7 deg. 

One way to solve the problems arising from the fact that 
| AC l | is small, | AC L | = 0.26, is to fix the angle of attack at 
a = 17 deg, corresponding to = -0.38, and use the bank angle to 
control both the longitudinal motion and the lateral motion 
(Refs. 8-9). Let L p = Lcosy denote the vertical projection of the lift, 
and let C LP = C L cosy denote the projected lift coefficient. Because 
— 1 < cosy < + 1, it is easy to see that = ^.76. Hence, the 

projected lift coefficient range is three times the lift coefficient 
range. To sum up, the control of the longitudinal motion is obtained 
by changing the modulus of the bank angle, while the control of the 
lateral motion is obtained by changing the sign of the bank angle. 
Reference 8 provides a nominal trajectory in which the bank angle 
history is represented by five constant-bank-angle segments and 
therefore involves four switches of the bank angle. 

1.2. AFE Optimal Trajectories . Independently of the present 
plans for the aeroassisted flight experiment, it is of interest 
to determine the optimal trajectories of the AFE spacecraft. This 
study was carried out in Refs. 10-14 under the following 
assumptions: (a) the angle of attack is constant, a = 17 deg; 
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(b) the spacecraft is controlled via the angle of bank; (c) the 
entry conditions are identical with those of the nominal 
trajectory of Ref. 8. In particular, in an inertial reference 

r 

frame, the entry path inclination is fixed, Yg = -4.49 deg. 

From the extensive numerical computations, the following conclusions 
were obtained in Refs. 10-14: 

(i) the optimal trajectories are two-subarc trajectories, 
with the bank angle constant in each subarc; hence, the control 
is bang-bang; 

(ii) in the atmospheric entry phase, the bank angle is 

y = 176.7 deg, yielding a positive vertical component of the lift, 
which in turn helps the path inclination to increase gradually from 
the entry negative value to nearly zero value; 

(iii) in the atmospheric exit phase, the bank anale is 

y = 5.5 deg, yielding a negative vertical component of the lift, which 
offsets the centrifugal force effects due to the curvature of the Earth, 
so as to ensure exit conditions compatible with the specified LEO; 

(iv) the horizontal component of the lift during the atmospheric 
entry phase and the horizontal component of the lift during the 
atmospheric exit phase have the same sign and the same order of 
magnitude; they are directed in such a way that they nearly offset 
the effects due to the Earth's rotation; in this way, the 
instantaneous orbital plane is almost identical with the initial 
orbital plane, meaning that the wedge angle is nearly zero during 

the atmospheric pass; this means that, for efficient flight, the 
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motion of the AFE spacecraft is nearly planar in an inertial 
space; in other words, one must avoid energy dissipation 
associated with the lateral motion. 

1.3. Present Research . This paper extends the optimization 
studies initiated in Ref. 10-14 to the case where (a) the AFE 
spacecraft is controlled via both the angle of attack and 
the angle of bank and (b) the entry conditions are identical 
with those of the nominal trajectory of Ref. 8, except the 
path inclination and the velocity, which are optimized. Optimal 
trajectories are studied for two transfer maneuvers: (DA) direct 

ascent to LEO; (IA) indirect ascent to LEO via PEO. 

While the motion of the AFE spacecraft in a 3D-space is 
described by a system of six ODEs, substantial simplifications 
are possible if one exploits these facts: (i) the instantaneous 

orbital plane is nearly identical with the initial orbital plane; 
(ii) the bank angle is small; and (iii) the Earth's angular 
velocity is relatively small. Under these assumptions, the 
complete system can be decoupled into two subsystems, one 
describing the longitudinal motion and one describing the lateral 
motion . 

The angle of attack history, the entry path inclination, and 
the flight time are determined via the longitudinal motion 
subsystem; in this subsystem, the total characteristic velocity 
is minimized subject to the specified LEO requirement. The angle 
of bank history is determined via the lateral motion subsystem; 
in this subsystem, the difference between the instantaneous bank 
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angle and a constant bank angle is minimized in the least 
square sense subject to the specified orbital inclination 
requirement . 

It is shown that both the angle of attack and the angle of 
bank are constant. This result has considerable importance in the 
design of nominal trajectories to be used in the guidance of AFE 
and AOT vehicles. 

1.4. Outline . Section 2 contains the system description, and 
Section 3 describes the optimal control problems to be solved. 
Section 4 introduces the decomposition technique, leading to the 
decoupling of the longitudinal motion (Section 5) from the lateral 
motion (Section 6). The experimental data are given in Section 7, 
and the numerical results are shown in Section 8. Finally, the 
conclusions are given in Section 9. 
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2 . System Description 

The motion of the AFE spacecraft takes place partly in 
space and partly in the atmosphere. For the purposes of this 
report, the trajectory begins at GEO (h = 19323 NM = 35786 km) 
and ends at LEO (h = 178 NM = 330 km) . It includes a preatmospheric 
branch, an atmospheric branch, and a postatmospheric branch. 
Depending on the nature of the postatmospheric branch, we consider 
two transfer maneuvers: (DA) direct ascent to LEO; (IA) indirect 
ascent to LEO via an intermediate parking Earth orbit (PEO, 
h = 197 NM = 365 km) . We assume that GEO, LEO, and PEO are 
circular orbits. 

For Transfer (DA), the key points of the maneuver are these: 
point 00, exit from GEO; point 0, atmospheric entry; point 1, 
atmospheric exit; point 11, entry into LEO. Point 00 is the 
apogee of the preatmospheric transfer orbit 00 -* 0; point 11 
is the apogee of the postatmospheric transfer orbit 1 ->• 11. 
Propulsive impulses are applied at two points: at point 00 
to deorbit from GEO; at point 11 in order to circularize the 
motion into LEO. See Fig. 2A. 

For Transfer (IA) , the key points of the maneuver are these: 
point 00, exit from GEO; point 0, atmospheric entry; point 1, 
atmospheric exit; point 22, entry into PEO; and point 11, entry 
into LEO. Point 00 is the apogee of the preatmospheric transfer 
orbit 00 0; point 22 is the apogee of both the first post- 

atmospheric transfer orbit 1 •> 22 and the second postatmospheric 
transfer orbit 22 -> 11; point 11 is the perigee of the second 
postatmospheric transfer orbit 22 -*• 11. Propulsive impulses are 
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applied at three points: at point 00 to deorbit from GEO; at 
point 22 in order to raise the height of the perigee of the second 

postatmospheric transfer orbit; and at point 11 in order to 
circularize the motion into LEO. See Fig. 2B. 

For the atmospheric portion (h < h ) of the trajectory of 

— a 

the AFE spacecraft, we employ an Earth-fixed system; for the 
space portion of the trajectory (h > h^) , we employ an inertial 
system; here, h = 400000 ft = 122 km denotes the thickness of 

cl 

the atmosphere. For h £ h a , we compute the air density using 
the US Standard Atmosphere, 1976 (Ref. 15); for h > h^ , we 
assume that the air density is zero. For both the atmospheric 
portion and the space portions of the trajectory, we neglect 
the effects due to the oblateness of the Earth; we assume that 
the gravitational field is central and obeys the inverse square 
law. 

2.1. Atmospheric Pass . With reference to the atmospheric 
portion of the trajectory of the AFE vehicle, the following 
additional hypotheses are employed: (a) the atmospheric pass is 

made with engine shut off; hence, the AFE spacecraft behaves as 
a particle of constant mass; (b) under extreme hypersonic 
conditions, the dependence of the aerodynamic coefficients on 
the Mach number and the Reynolds number is disregarded; (c) the 
sideslip angle is zero; hence, the side force component of the 
aerodynamic force is zero; (d) the AFE spacecraft is controlled 
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via the angle of attack and the angle of bank. 

2.2. Differential System . With the above assumptions, 
upon using an Earth-fixed system, and upon normalizing the 
interval of integration to unity, the equations of motion include 
the kinematical equations (Ref. 10) 

9 = x (Vcosy cosy/rcos<{>) , (l a ) 

<p = i (-Vcosy siny/r) , (lb) 

r = x (Vsiny) , (1°) 

and the dynamical equations (Ref. 10) 

V = x (-D/m - gsiny) 

2 2 

4- x [ m r(siny cos + cosy siny cos<J) simji) ] , (2a) 

y = x [ (L/mV) cosy + (V/r - g/V)cosy + 2m cosy cos<}>] 

2 2 

+ x [ (to r/V) (cosy cos <)> - siny siny cost}) sin<J>)] f (2b) 

X = x [ (L/mV) siny/cosy + (V/r) cosy cosy tanij>] 

+ x [2to (sini}> + tany siny cost}))] 

2 

+ x [(to r/V) cosy cost}) sin<}>/cosy] . (2c) 

In the dynamical equations, the symbol w denotes the angular 
velocity of the Earth; terms linear in m are due to the Coriolis 
acceleration; terms quadratic in m are due to the transport 
acceleration. Also in the dynamical equations, the local 
acceleration of gravity is given by 
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g = y e /r 2 , (3) 

where y g denotes the Earth's gravitational constant. In addition, 

* the aerodynamic forces are given by 

D = (l/2)C D (ot) p (h)SV 2 , (4a) 

L = (l/2)C L (a)p(h)SV 2 , (4b) 

where the air density p depends on the altitude h, with 

h = r - r . (4c) 

e 

For given initial conditions, parameter t, and controls a(t) 
and y(t), Eqs. ( 1 ) - ( 4 ) can be integrated in forward time over the 
time interval 0 <_ t <_ 1. Here, the initial time t = 0 corresponds 
to atmospheric entry and the final time t = 1 corresponds to 
atmospheric exit. 

2.3. Control Constraint . To obtain realistic solutions, 
the presence of upper and lower bounds on the angle of attack 
is necessary. Therefore, the two-sided inequality constraint 

a < a < a rT , 0 £ t £ 1, (5) 

L u 

must be satisfied everywhere along the interval of integration. 

2.4. Transformation Relations . General transformation 
relations allow one to pass from quantities computed in an 
Earth-fixed system to quantities computed in an inertial system 
(direct relations) , and viceversa (inverse relations). See Refs. 10-11 


for details. 
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For the kinematic variables, the direct relations are 
6 = 6 + Toot, (6a) 

<p = <P, (6b) 

r = r, (6c) 

and the inverse relations are 

0 = 0- Toot, (7a) 

<j> = cj>, (7b) 

r = r. (7c) 

For the dynamic variables, the general transformation 
relations are 

Vcosy cosy = Vcosy cosy + wr cos4>, (8a) 

Vcosy siny = Vcosy siny, (8b) 

Vsiny = Vsiny, (8c) 

and they imply the direct relations 

2 2 
V = /[V + 2oorVcosy cosy cos4> + (cor cosc(>) ] , (9a) 

2 2 
tany = Vsiny/ /[ (Vcosy) + 2corVcosy cosy cos<j) + (cor cos4>) ], (9b) 

tany = Vcosy siny/(Vcosy cosy + oorcos<|>) , (9c) 


and the inverse relations 
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~ 2 — ~ ~ - 2 

V = /[V - 2cjrVcosy cosx cost)) + (oorcosi) ], (10a) 

tany = Vsiny//[ (Vcosy) - 2wrVcosy cosy cosij) + (air cost})) ] , (10b) 

tany = Vcosy sinx/(Vcosy cosy - wrcostji) . (10c) 

If the Earth's rotation is neglected (w = 0), the inertial 
quantities appearing in Eqs.(6)~(10) become identical to the 
corresponding Earth-fixed quantities. However., this is not the 
case if the Earth's rotation is considered (co ^ 0) . 

2.5. Orbital Elements . Once the state variables are 
known, one can compute some important quantities such as the 
orbital inclination, the longitude of the ascending node, and 
the wedge angle; the latter is the angle between the instantaneous 
orbital plane and the entry orbital plane. In the inertial 
system, these quantities are defined through the relations 

cosi = cos<j> cosy, (11a) 

sin(0 - Q) = coti tan<}> , (lib) 

cosn = sini sinig cos(0 - fig) + cosi cosig. (11c) 

In the Earth -fixed system, the relations analogous to (11) are 
defined below: 

cosi = cosij) cosy, (12a) 

sin(Q - (7) = coti tant}», (12b) 

cosg = sini sinig cos(ft - fig) + cosi cosig. 


(12c) 
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If the Earth's rotation is neglected (to = 0) , the inertial 
quantities appearing in Eqs. (11) become identical to the 
corresponding Earth-fixed quantities. However, this is not the 
case if the Earth's rotation is considered (w ^ 0). 

2.6. Initial Conditions . At atmospheric entry, the initial 
values of the state variables 0 Q , <>q, r Q , x 0 are gi ven in 

the inertial system. In particular, 

r _ r . (13a) 

0 a 

Because of the relations (11) , the initial values of the orbital 
elements i^, Oq are known. In particular, 

0 Q = 0, d3b) 

by definition. The initial values of the state variables V Q , y Q 
must be consistent with the relation 

r 00 (2v " - 5 0> - 2r 00 r a V * + r a^0 cos2 ^0 = °' (13c) 

where V = V = /7y /r )= 7.831 km/sec is a reference velocity, i.e., 
* B. 0 3 

the circular velocity at r = r . This relation arises from 
energy conservation and angular momentum conservation applied 
to the preatmospheric transfer orbit 00 -*■ 0 connecting GEO with 
atmospheric entry. 

2.7. Final Conditions . At atmospheric exit, the final 
time t is free and must be determined in such a way that 

r = r (14a) 

r l a 

The final values of the state variables 0^, 4^, X]_ must be 
consistent with the relation 
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n 1 = o. (14b) 

This means that the final value of the wedge angle must vanish, 
hence that the exit orbital plane must be identical with the 
entry orbital plane. The final values of the state variables 
V l' Y 1 raust be cons i stent with either the relation 

(DA) r^(2vj - V 2 ) - 2r i;i r _V 2 + r^cos 2 y 1 = 0 (14c) 

or the relation 

(IA) r^ 2 (2vf - vf) - ? r 22 r a V * + r a* 1 cos ^l = °‘ (14d) 

Equation (14c) refers to the direct ascent case and arises 
from energy conservation and angular momentum conservation 
applied to the postatmospheric transfer orbit 1 ■* 11 connecting 
atmospheric exit with LEO. Equation (14d) refers to the indirect 
ascent case and arises from energy conservation and angular 
momentum conservation applied to the postatmospheric transfer 
orbit 1 ^ 22 connecting atmospheric exit with PEO. 

2.8. Summary . The relations governing the atmospheric pass 
include the differential system (l)-(4), the control constraint 
(5), and the boundary conditions (13)-(14). In the boundary 
conditions, the inertial quantities are related to the Earth- 
fixed quantities via the transformation relations (6) -(10); 
also, the inertial quantities are related to the orbital elements 
via the trigonometric relations (11) . In this formulation, 
the independent variable is the time t, 0 <_ t <_ 1. The dependent 
variables include six state variables [0(t), <}>(t), r(t), V(t), 

Y ( t) , two control variables [a(t), y (t) ] , and one 

parameter (t) . 
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3 . Optimal Control Problems 

Subject to the previous constraints, different optimal 
control problems can be formulated, depending on the performance 
index chosen and the type of transfer maneuver considered. Here, 
we focus attention on the minimization of the total characteristic 
velocity AV, which is a measure of the propellant 
required for orbital transfer. 

Problem (DA ) . This problem refers to the direct ascent 
to LEO. The functional to be minimized is given by 

I = AV = AV Q0 + AV U , (15a) 

with 


4V 00 = /(r a /r 00 )V * - (r a /r 00 )V 0 COSY 0' 


(15b) 


AV 11 = / < r a /r ll ,V * - (r a /r ll )V l COSY l' 


(15c) 


Problem (IA ) . This problem refers to the indirect ascent 
to LEO via PEO. The functional to be minimized is given by 


I - 4V - 4V 00 + AV 22 + 4V 11’ 


(16a) 


with 


A ^00 = / ’ (r a /r 00 ,V * " (r a /r 00 } V OS V 


(16b) 


AV 22 = ^[2r a r 11 /(r 11 r 22 + r* 2 ) ]V # -(^/r^JV^os^, (16c) 


A ^ll = /l2r a r 22 /(r il r 22 + r il )]V * " / ‘ (r a /r ll )V *- 


(16d) 
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4 . Decomposition Technique 

The study of the optimal trajectories of an AFE vehicle 
can be simplified to a considerable degree if the complete 
system of differential equations, inequality constraints, and 
boundary conditions can be decomposed into a subsystem governing 
the longitudinal motion [state variables r(t), V(t), y(t); 
control variable a(t); unknown parameter x] and a subsystem 
governing the lateral motion [state variables 0(t), x(t); 

control variable y(t); known parameter x] . The decomposition 
is possible if the following approximations are introduced into 
the equations of motion: 


1 = V 

(17a) 

y 2 << l. 

(17b) 

CO 2 = 0. 

(17c) 


Approximation (17a) means that the instantaneous orbital plane 
is nearly identical with the initial orbital plane. Approximation 
(17b) means that the bank angle is small. Approximation (17c) 
means that terms quadratic in to are small with respect to terms 
linear in w and terms not containing to. The following relations 
arise as a consequence of the approximations (17): 


cosig = cose}) cosx, 

(18a) 

cosy = 1, 

(18b) 
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~ 2 

V = /”[V + 2wrVcosy cosig] , 

~ ~ ~ 2 2 

Vcosy = /[V cos y + 2u)rVcosy cosig] . 

Relation (18a) is due to (12a) and (17a) . Relation (18b) is 
a restatement of (17b) . Relation (18c) is obtained by combining 
(9a), (17c), and (18a). Relation (18d) arises from (8a), (8b), 

(17c), and (18a). As the subsequent analysis shows, the 
approximations (17) and the implications (18) are essential 
to a simplified study of the longitudinal motion (Section 5) , 
but nonessential, hence optional, to a simplified study of 
the lateral motion (Section 6) . 


(18c) 

(18d) 
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5 . Longitudinal Motion 

In the Earth-fixed system, the subsystem governing the 
longitudinal motion includes the differential equations 

r = x [Vsiny] , (19a) 

V = x [-D/m - gsiny] , (19b) 

y = x [L/mV + (V/r - g/V)cosy + 2wcosi Q ] , (19c) 

the control inequality constraint 


a L i a < V 


( 20 ) 


the initial conditions 


r o = V 


(21a) 


2V * (r 00 " r 00 r a ) + 2 “ r a V 0 COSY 0 cosi 0 (r a 


2 x 

r oo } 


+ V^(r^cos 2 Y 0 - r 2 0 ) = 0, 


(21b) 


the final conditions 


r, = r 
1 a 


(22a) 


and either 


(DA) 2vl(r 2 ±1 - r n r a ) + 2 t or.V 1 cosY 1 cosier* - r^) 


+ V 2 (r 2 cos 2 y 1 - r^) = 0, 


(22b) 
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or 

(IA) 2V*(r 2 22 - r 22 r a ) + 2u3r a V ] cosY 1 cosier* - r^) 

+ V^(r a cos^Y-^ - ^22^ = (22c) 

Relations (19) -(22) are obtained from (lc) , (2a), (2b), (5), 

(13a), (13c), (14a), (14c), (14d) after accounting for (6), (8), the 

approximations (17), and the implications (18). 

5.1. Performance Indexes . Using (17)-(18), the performance 
indexes of the optimal control problems of Section 3 can be 
reformulated below. For the direct ascent to LEO, Eqs . (15) become 

(DA) I = AV = AV q0 + AV U , < 23a) 

with 


w 00 = / <V r 00 )V * 

- < r a /r 00 ,/7v 0 oos2Y 0 + 2 “ r a V 0 COSY 0 cosi 0>' 
AVu = /(r a /r 11 )V* 

- 2 2 

- (r a /r i;L ) /(V 1 cos Y x + 2wr a V l COSY l cosi Q ) . 


For the indirect ascent to LEO via PEO, Eqs. (16) become 


(IA) I = AV = AV Q0 + AV 22 + AV i;l , 
with 

A ^00 = ^ r a /r 00 )V * 

* (r a /r 00 )/ < V 0 cos2Y 0 + 2 “ r a V 0 COSY 0 cosl 0 ) ' 


(23b) 


(23c) 


(24a) 


(24b) 
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AV 


22 


/ U r a r ll /(r U r 22 + r 22 )]v, 


9 9 

(r a /r 22 ) /(V£cos y^^ + 2wr a V 1 cosy 1 cosi Q ) 


AV 


11 


v/ '- 2r a r 22 // ^ r ll r 22 + r ll^ V *~ v/ ^ r a //r ll^ V * 


(24c) 
( 24d) 


5.2. Summary . The relations governing the 
longitudinal motion include the differential system (19), the 
control constraint (20), and the boundary conditions (21) -(22). 

In this formulation, the independent variable is the time t, 

0 < t < 1. The dependent variables include three state variables 
[r(t), V(t), y (t) ] , one control variable [a(t)], and one 
parameter (t) . These variables must be determined in such a way 
that the performance index (23) is minimized in the direct ascent 
case and the performance index (24) is minimized in the indirect 


ascent case. 
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6 . Lateral Motion 

In the Earth-fixed system, the subsystem governing the 

lateral motion includes the differential equations 




• 

0 = x (Vcosy cosy/rcos4>) , 

(25a) 

— 


« 

<j> = t (-Vcosy sinx/r) , 

(25b) 

— 


X = x [ (L/mV) siny/cosy + (V/r)cosy cosy tan<}>] 


- - 


+ T[2u>(sin<{> + tany siny cost}))] 




2 

+ x [(w r/V) cosy cos<|> sin^/cosy] , 

(25c) 


the 

initial conditions 




e Q = given, 

(26a) 



4> 0 = given. 

(26b) 



X 0 = given, 

(26c) 

— 

and 

the final condition 


— 


r\ 1 = 0. 

(27a) 

— 

are 

It must be noted that the functions r (t) , V(t), y (t) , a(t), x 
known from the solution of the optimization problem associated 


with the longitudinal motion (Section 5) . It must also be noted 
that, on account of the definitions (11) and the transformation 
relations (6) and (9) , the final condition (27a) can be 
rewritten in the functional form 
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Tl-L ( ©-L / 4> 1 * X x ) = 0, (27b) 

which replaces (27a) . 

6.1. Performance Index . In the subsystem (25) -(27), the 
dependent variables include three state variables [9(t), <t>(t), 
x(t)] and one control variable [y(t)]. There is no parameter, 
since x is known. Therefore, the subsystem (25)- (27) admits an 
infinite number of solutions for the bank angle y (t) . The 
solution can be rendered unique by requiring that an optimization 
criterion be met, such as minimizing the functional 


J 0 

where n is a parameter. The solution of the optimal control 
problem (25) -(28) has the form 


y (t) = tt = const . 


(29) 


Indeed, it can be readily shown that the multipliers associated 
with the constraints (25) -(27) all vanish. 

6.2. Two-Point Boundary-Value Problem . For the purposes 
of explanation, observe that (29) implies that 


y = 0 . 


(30) 


Then, the subsystem (25) -(27) augmented by (30) constitutes a 
two-point boundary-value problem (TPBVP) in which the unknowns 
are the functions 0 (t) , <p (t) , x(t), y(t). Since there are four 
differential equations and four boundary conditions, one surmises 
that a solution might exist; however, because the TPBVP (25) -(27) 
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and (30) is nonlinear, the existence of a solution must be 

confirmed by numerical tests. 

6.3. Remark . While the hypothesis (17c) is essential 

to a simplified study of the longitudinal motion, it is optional 

to a simplified study of the lateral motion. Indeed, the 

previous considerations stand unchanged regardless of whether 
2 

the to term is dropped or retained in (25c) . If it is 
retained, then a slightly more precise solution is obtained 
for the lateral motion. 
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7 • Experimental Data 

The following data were used in the numerical experiments. 

Spacecraft Data . The AFE configuration is shown in Fig. 1. 

The mass of the AFE spacecraft is m = 1678 kg; the reference 

surface area is S = 14.31 m . For this configuration, the drag 

coefficient = C^(a) is shown in Fig. 3A, the lift coefficient 

= < “ L ( a ) i s shown in Fig. 3B, the lift-to-drag ratio modulus 

E = E(cc) is shown in Fig. 3C, and the drag polar C = C (C ) 

D D L 

is shown in Fig. 3D. See also Table 1. 

The functions appearing in Fig. 3 can be approximated 
by the relations 


C D = A Q + A-^a + A 2 a , 


C L = B 0 + B^a + B 2 a ' 


(31a) 

(31b) 


with 


E * 

C Q + C,a 

+C 2 a^ , 




(31c) 

C D = 

D 0 + D 1 C 

L + °2 C L 2 

/ 



Old) 

> 

o 

II 

1.565, 

A 1 = ' 

-0.306 , 

A 2 

= -1.375, 

(32a) 

B o = 

-0.036, 

B 1 = ' 

-1.557, 

B 2 

= 1.344, 

(32b) 

II 

o 

U 

0.036, 

o 

II 

0.836, 

C 2 

= -0.041, 

(32c) 

II 

o 

Q 

1.235, 

D l 

-2.379, 

D 2 

= -5.473. 

(32d) 
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The angle of attack is subject to Ineq. (5) , 


a L - a - V 


corresponding to 


C LL - C L - C LU" 


(33a) 


(33b) 


Note that the lower bound in (33a) corresponds to the upper 
bound in (33b), and viceversa. This is due to the fact that, 
because of the raked-cone configuration of the AFE spacecraft, 
the lift coefficient is a decreasing function of the angle of 
attack (see Fig. 3B) . 

In the numerical experiments, three alternative angle of 
attack ranges were considered, hence three alternative lift 
coefficient ranges were considered , through the following choices: 


(i) 

a L = 

7.0 deg. 

a u = 

27.0 deg. 

(34a) 

(ii) 

a L = 

7.0 deg, 

a u = 

17.0 deg. 

(34b) 

(iii) 

a L = 

7.0 deg, 

a u = 

10.0 deg. 

(34c) 

corresponding to 





(i) 

C LL 

r* 

• 

o 

1 

II 

C LU 

= -0.21, 

(35a) 

(ii) 

C LL 

II 

1 

o 

• 

CO 

C LU 

= -0.21, 

(35b) 

(iii) 

C LL 

= -0.27, 

C LU 

= -0.21. 

(35c) 


For these cases, the values of the maximum lift-to-drag ratio 


\ 


modulus are 
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(i) 

E 

= 0.42, 

max 


(ii) 

E 

11 

o 

ro 

00 

max 


(iii) 

E 

max 

II 

o 

M 

00 


(36a) 

(36b) 

(36c) 


and occur for C = C . 

Li LL 

Physical Constants . The radius of the Earth is r g = 6378 km; 
the radius of the outer edge of the atmosphere is r = 6500 km; 

cl 

the thickness of the atmosphere is h^ = 122 km; the Earth's 

3 2 

gravitational constant is = 0.3986E+06 km /sec ; the circular 
velocity at r = r is V = 7.831 km/sec; the angular velocity 

3. a 

of the Earth is w = 0.7292E-04 rad/sec. 

Atmospheric Model . The assumed atmospheric model is that of 
the US Standard Atmosphere, 1976 (Ref. 15). In this model, the 
values of the density are tabulated at discrete altitudes. For 
intermediate altitudes, the density is computed by assuming an 
exponential fit for the function p (h) . 

Dynamic Pressure . The dynamic pressure is computed with 
the formula 


DP = (1/2) pV 2 . 


(37a) 


Heating Rate . The stagnation point heating rate is computed 
with the formula 


3.07 


HR = c/(p/P*) (V/V*) 


(37b) 
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3 

Here, p* = 0.3097E-03 kg/m is a reference density, the density 
at h. = 60 km; V. = V = 7.831 km/sec is a reference velocity; 

x * 3. 

2 

the constant C = 282.3 W/cm represents the stagnation point 
heating rate at p = p* and V = V*, based on a nose radius of 
one foot. 

Atmospheric Entry . In the inertial system, the given 
initial conditions are as follows: the longitude is 0 q = 

-134.52 deg; the latitude is = -4.49 deg; the altitude is 
h Q = 122 km, corresponding to the radius r Q = 6500 km; the 

heading angle is Xg = ~28.13 deg. The orbital inclination is 
i = 28.45 deg; the longitude of the ascending node is = 

-126.19 deg; the wedge angle is = 0.00 deg. 

Atmospheric Exit . In the inertial system, the desired 
final conditions are as follows: the altitude is h^ = 122 km, 
corresponding to the radius r 1 = 6500 km; the orbital inclination 
is i = 28.45 deg; the longitude of the ascending node is 
= -126.19 deg; the wedge angle is = 0.00 deg. 

Transfer (DA ) . This is the direct ascent to LEO. The GEO 
conditions are as follows: the altitude is h Q0 = 35786 km, 
corresponding to the radius r^ = 42164 km; the path inclination 

is y qo = 0.00 deg. The LEO conditions are as follows: the 
altitude is h = 330 km, corresponding to the radius r ^ = 

6708 km; the path inclination is =0.00 deg. 


\ 
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Transfer (IA) . This is the indirect ascent to LEO via PEO. 

The GEO conditions are as follows: the altitude is hgg = 35786 km, 
corresponding to the radius r^g = 42164 km; the path inclination 

is Ygg = 0.00 deg. The PEO conditions are as follows: the 
altitude is = 365 km, corresponding to the radius r ^ = 

6743 km; the path inclination is y ^ = 0-00 deg. The LEO 
conditions are as follows: the altitude is h^ = 330 km, 

corresponding to the radius r^ = 6708 km; the path inclination 
is y^i = 0.00 deg. 


» 
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8 . Numerical Results 

The decomposition technique of Sections 4-6 was implemented 
for the AFE spacecraft in connection with the experimental data 
of Section 7. First, the longitudinal motion subsystem was 
solved with respect to the functions r(t), V(t) , y(t), a(t) 
and the parameter t which minimize the total characteristic 
velocity, while satisfying the desired LEO requirement; here, 
the sequential gradient-restoration algorithm (SGRA, see Refs. 16-18 
for recent versions) was employed. Then, with the above functions 
and parameter fixed, the lateral motion subsystem was solved 
with respect to the functions 0 (t) , <f>(t), x(t)» v(t) which 
minimize the difference y(t) - it in the least square sense, tt 
being a parameter. This leads to a nonlinear two-point boundary- 
value problem (TPBVP) ; here, the modified quasilinearization 
algorithm was employed in conjunction with the method of particular 
solutions (MQA/MPS, see Refs. 19-21). 

8.1. Effect of the Control Bounds . First, (DA) the direct 
ascent to LEO was considered, and optimal trajectories (OT) 
were computed for three alternative angle of attack ranges, 
corresponding to three alternative lift coefficient ranges, 
described by (33) -(36). In each case, consistently with the 
results of Ref. 7, it was found that optimality in the longitudinal 
motion is achieved by flying at the lift coefficient lower bound, 
hence at the angle of attack upper bound. Concerning the 
lateral motion, the least square deviation of the bank angle from 
a constant value is achieved by flying at constant bank angle. 
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Summary results can be found in Table 2 , which shows the 
following quantities for DAOT Cases (i) , (ii) , (iii) : the number 
of bank angle switches, the flight time, the angle of attack, the 
lift coefficient, and the angle of bank; the minimum altitude, 
the peak dynamic pressure, and the peak heating rate; the 
entry path inclination, the peak change of orbital inclination, 
the peak change of longitude of the ascending node, and the 
peak wedge angle; the characteristic velocity components and 
the total characteristic velocity. 

As can be seen, the total characteristic velocities of 
DAOT Cases (i) , (ii) , (iii) are nearly the same. With respect 
to the characteristic velocity of Case (ii) , that of Case (i) 
is 0.1% lower, while that of Case (iii) is 0.2% higher. However, 
the solution (ii) is to be preferred because of the following 
consideration: it is characterized by C M = 0 , where C^ is the 
moment coefficient (Ref. 8), while this is not the case with 
solutions (i) and (iii). 

Comparison of DAOT Case (ii) with the DAOT solution presented 
in Ref. 14 shows that the entry path inclination is now flatter 
(by 0.36 deg). Consequently, the minimum altitude is higher (by 
about 0.6 km); hence, the peak dynamic pressure and the peak 
heating rate are lower than the corresponding quantities in Ref. 14. 

8.2. Comparison of AFE Trajectories . Based on the results 
of the previous section, we restrict the analysis to Case (ii) , 
namely the case where Ineqs. (33) have the form 
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(ii) 

7.0 < a < 17.0 deg , 

(38a) 

(ii) 

-0.38 < C < -0.21. 

— Ij — 

(38b) 


Both (DA) the direct ascent to LEO and (IA) the indirect 
ascent to LEO via PEO were considered, and optimal trajectories 
(OT) were computed. Once more, it was found that optimality 
in the longitudinal motion is achieved by flying at the upper 
bound angle of attack, hence at the lower bound lift coefficient. 
Concerning the lateral motion, optimality is achieved by flying 
at constant bank angle. 

For comparison purposes and for (IA) the indirect ascent to 
LEO via PEO, a reference trajectory (RT) was also computed. 

This is the nominal trajectory of Ref. 8, which is flown at 
a = 17.0 deg using N g = 4 switches of the bank angle, hence 
N = 5 constant values of the bank angle, specif ically : 

]i ± = -180.0, v 2 = -43.9, y 3 = 89.9, y 4 = -89.9, y 5 = 0.0 deg. (39) 

For trajectories DAOT , IAOT, and I ART , Table 3 shows the 
following quantities: the number of bank angle switches, the 
flight time, the angle of attack, the lift coefficient, and 
the angle of bank; the minimum altitude, the peak dynamic pressure, 
and the peak heating rate; the entry path inclination, the peak 
change of orbital inclination, the peak change of longitude of 
the ascending node, and the peak wedge angle; the characteristic 
velocity components and the total characteristic velocity. 
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For trajectories DAOT, IAOT, and IART, and for both the 
Earth-fixed system and the inertial system. Tables 4-5 show 
the entry values (Table 4) and the exit values (Table 5) 
of the following quantities: longitude, latitude, altitude, 
radius, velocity, path inclination, and heading angle; orbital 
inclination, longitude of the ascending node, and wedge angle. 

For trajectories DAOT, IAOT, and IART, Fig. 4 shows the 
time histories of the following quantities: relative longitude 
(Fig. 4A) , latitude (Fig. 4B) , altitude (Fig. 4C) , relative 
velocity (Fig. 4D) , relative path inclination (Fig. 4E) , 
relative heading angle (Fig. 4F) ; angle of attack (Fig. 4G) , 
lift coefficient (Fig. 4H) , bank angle (Fig. 41) ; dynamic pressure 
(Fig. 4J) , heating rate (Fig. 4K) ; inertial orbital inclination 
(Fig. 4L) , inertial longitude of the ascending node (Fig. 4M) , 
and inertial wedge angle (Fig. 4N) . 

From Tables 3-5 and Fig. 4, the following comments arise: 

(a) The total characteristic velocity of trajectory DAOT 
is less than that of trajectory IAOT, which in turn is less 
than that of trajectory IART. 

(b) Compared with trajectory IART, trajectories DAOT 

and IAOT are characterized by flatter entry (by 0.36 deg) , hence 
higher minimum altitude (by 3.3 km), lower peak dynamic pressure, 
lower peak heating rate, and longer flight time. 

(c) The longer flight time of trajectories DAOT and IAOT 
vis-a-vis trajectory IART can be explained as follows: while 
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the velocity depletion V Q - is about the same for the 
three trajectories, the first two are flown at higher altitude, 
hence smaller air density, hence less aerodynamic drag. 

(d) Compared with trajectory IART, trajectories DAOT 
and IAOT involve much smaller values of the peak change of 
orbital inclination, peak change of longitude of the ascending 
node, and peak wedge angle. Indeed, for trajectories DAOT and 
IAOT, these quantities are nearly zero, within 1/100 deg. On 
the other hand, for trajectory IART, max|Ai| = 1.08deg, 
max | Aft | = 0.31 deg, max(n) = 1.08 deg. 

8.3. Control Considerations . Trajectories DAOT, IAOT, 
and IART are all flown at the same angle of attack a = 17.0 deg, 
hence at the same lift coefficient C^ = -0.38. This being the 
case, the control margin for the angle of attack, hence the 
control margin for the lift coefficient, is the same for the 
three trajectories. 

Concerning the angle of bank, trajectories DAOT and IAOT 
are flown with y = 5.14 deg, while trajectory IART is flown 
with five subsequent values of the bank angle [see (39)]. Hence, 
trajectories DAOT and IAOT have a larger control margin for 
the angle of bank than trajectory IART. 

It should be pointed out that trajectories DAOT, IAOT, 
and IART are inherently unstable if flown open-loop. For 
example, suppose that, at time instant t, there is a positive 
difference between the flight velocity and the nominal velocity. 
The higher flight velocity yields a higher centrifugal force 
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due to the Earth's curvature, hence a higher value of y, hence a 
higher altitude, hence a lower density, hence a lower drag, hence 
less velocity decrease vis-a-vis the nominal trajectory. 
Therefore, at time instant t + At, the positive difference 
between the flight velocity and the nominal trajectory has 
further increased, and so on. To sum up, (i) guidance and 
control systems must be designed with care; and (ii) feedback 
control schemes must be developed to ensure the stability of 
the AFE trajectory. 
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9 . Conclusions 

In this report, the GEO-to-LEO transfer of an AFE spacecraft is 
considered, and optimal trajectories are determined by minimizing 
the total characteristic velocity. The optimization is performed 
with respect to the time history of the controls (angle of attack 
and angle of bank) , the entry path inclination and the flight 
time being free. Two transfer maneuvers are considered: (DA) 

direct ascent to LEO; (IA) indirect ascent to LEO via PEO. 

While the motion of the AFE spacecraft in a 3D-space is 
described by a system of six ODEs, substantial simplifications 
are possible if one exploits the smallness of three key 
quantities: the change of orbital inclination; the bank angle; 
and the Earth's angular velocity. Indeed, the complete system can 
be decoupled into two subsystems, one describing the longitudinal 
motion and one describing the lateral motion. 

The angle of attack history, the entry path inclination, and 
the flight time are determined via the longitudinal motion 
subsystem; in this subsystem, the total characteristic velocity 
is minimized subject to the specified LEO requirement. The angle 
of bank history is determined via the lateral motion subsystem; 
in this subsystem, the difference between the instantaneous bank 
angle and a constant bank angle is minimized in the least square 
sense subject to the specified orbital inclination requirement. 

From the extensive numerical computations, the following 
conclusions were obtained: 
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(i) The optimal trajectories are one-subarc trajectories, 
with constant angle of attack and constant angle of bank. 
Specifically, a = 17.0 deg and y = 5.14 deg. 

(ii) Throughout the atmospheric pass, the (a, y) -pair 
yields a negative vertical component of the lift, which offsets 
the centrifugal force effects due to the curvature of the Earth, 
so as to ensure exit conditions compatible with the specified 
LEO. 

(iii) Throughout the atmospheric pass, the (a,y)-pair 
yields a horizontal component of the lift which nearly offsets 
the effects due to the Earth's rotation. In this way, the 
instantaneous orbital plane is almost identical with the initial 
orbital plane, meaning that the wedge angle is nearly zero. This 
means that, for efficient flight, the motion of the AFE 
spacecraft is nearly planar in an inertial space; in other words, 
one must avoid energy dissipation associated with the lateral 
motion. 

(iv) In an inertial reference frame, the entry path 
inclination is Yq = —4.13 deg, thus yielding trajectories about 
0 . 36 deg flatter than both the nominal trajectory of Ref. 8 and 
the optimal trajectories of Ref. 14. In turn, this results in 
higher minimum altitude, lower peak dynamic pressure, and lower 
peak heating rate. 
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Table 1. AFE spacecraft, aerodynamic data. 


— 

a (deg) 

C L 

C D 

C L /C D 

C M 

— 

7 

-0.2088 

1.5078 

-0.1385 

0.0295 


9 

-0.2473 

1.4826 

-0.1668 

0.0232 


11 

-0.2842 

1.4555 

-0.1953 

0.0171 


13 

-0.3187 

1.4258 

-0.2235 

0.0113 


15 

-0.3512 

1.3915 

-0.2524 

0.0054 

- 

17 

-0.3807 

1,3549 

-0.2810 

-0.0004 


19 

-0.4065 

1.3132 

-0.3095 

-0.0064 


21 

-0.4283 

1.2688 

-0.3376 

-0.0123 


23 

-0.4458 

1.2187 

-0.3658 

-0.0181 


25 

-0.4595 

1.1683 

-0.3933 

-0.0235 

— 

27 

-0.4707 

1.1182 

-0.4210 

-0.0287 
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Table 2 . 

Effect of the control bounds, 
main properties. 


Quantity 

DAOT 

DAOT 

DAOT 

Units 


Case (i) 

Case (ii) 

Case (iii) 


N s 

0 

0 

0 

- 

X 

1045 

845 

689 

sec 

a 

27.0 

17.0 

10.0 

deg 

C L 

1 

O 

-J 

1 

O 

u> 

CO 

r- 

CM 

• 

O 

1 

- 

V 

3.45 

5.14 

3.18 

deg 

min (h) 

78.6 

78.1 

77.3 

km 

max (DP) 

991 

1040 

1160 

N/m 2 

max (HR) 

134 

135 

141 

W/cm 2 

*0 

-4.08 

-4.13 

-4.20 

deg 

max j Ai | 

0.01 

0.01 

0.01 

deg 

max | Aft | 

0.01 

0.01 

0.01 

deg 

max(n) 

0.01 

0.01 

0.01 

deg 

^00 

1.490 

1.490 

1.490 

km/sec 

AVn 

0.070 

0.072 

0.075 

km/sec 

AV 

1.560 

1.562 

1.565 

km/ sec 
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Table 3. Comparison of AFE trajectories, 
main properties. 


Quantity 

DAOT 

IAOT 

I ART 

Units 


Case (ii) 

Case (ii) 



N s 

0 

0 

4 

- 

X 

845 

821 

486 

sec 

a 

17.0 

17.0 

17.0 

deg 

C L 

-0.38 

-0.38 

-0.38 


V 

5.14 

5.14 

(*) 

deg 

min (h) 

78.1 

78.1 

74.8 

km 

max (DP) 

1040 

1039 

1578 

N/m 2 

max (HR) 

135 

135 

156 

W/crn 2 


-4.13 

-4.13 

-4.49 

deg 

max | Ai | 

0.01 

0.01 

1.08 

deg 

max | Aft | 

0.01 

0.01 

0.31 

deg 

max ( n) 

0.01 

0.01 

1.08 

deg 

A ^00 

1.490 

1.490 

1.491 

km/sec 

av 22 

0.000 

0.072 

0.088 

km/sec 

AVn 

0.072 

0.010 

0.010 

km/sec 

AV 

1.562 

1.572 

1.589 

km/sec 

(*) Values 

of y^ , y 2 , 

y 3' y 4 ' y 5 are 

given by Eqs. 

(39) . 
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Table 4. Comparison of AFE trajectories, 
entry conditions. 


Quantity 

DAOT 

IAOT 

I ART 

Units 


Case (ii) 

Case (ii) 



e o 

-134.52 

-134.52 

-134.52 

deg 

CD i 
O 

-134.52 

-134.52 

-134.52 

deg 

^0 

-4.49 

-4.49 

-4.49 

deg 

h o 

122 

122 

122 

km 

r o 

6500 

6500 

6500 

km 

v o 

9.906 

9.906 

9.895 

km/sec 

^0 

10.308 

10.308 

10.308 

km/sec 

Y o 

-4.296 

-4.296 

-4.675 

deg 

X 0 

-4.128 

-4.128 

-4.487 

deg 

x o 

-29.42 

-29.42 

-29.42 

deg 

Xo 

-28.13 

-28.13 

-28.13 

deg 

1 0 

29.73 

29.73 

29.73 

deg 

1 o 

28.45 

28.45 

28.45 

deg 

“o 

-126.62 

-126.62 

-126.62 

deg 

5o 

-126.19 

-126.19 

-126.19 

deg 

n o 

0.00 

0.00 

0.00 

deg 

^0 

0.00 

0.00 

0.00 

deg 





46 


AAR-250 


Table 5. Comparison of AFE trajectories, 
exit conditions. 


Quantity 

DAOT 

IAOT 

I ART 

Units 


Case (ii) 

Case (ii) 



9 1 

-78.90 

-80.50 

-103.31 

deg 

9 1 

-75.37 

- 77.07 

-101.28 

deg 

*! 

22.79 

22.28 

12.84 

deg 

h i 

122 

122 

122 

km 

r l 

6500 

6500 

6500 

km 

V 1 

7.481 

7.491 

7.462 

km/sec 

*1 

7.882 

7.892 

7.876 

km/sec 

Y 1 

0.803 

0.860 

1.358 

deg 

X 1 

0.762 

0.816 

1.287 

deg 

X 1 

-18.52 

-19.21 

-27.15 

deg 

X 1 

-17.52 

-18.17 

-25.61 

deg 

A 1 

29.05 

29.10 

29.82 

deg 

± 1 

28.46 

28.46 

28.46 

deg 

°1 

-128.04 

-127.92 

-126.74 

deg 

& 1 

-126.20 

-126.20 

-126.16 

deg 

n l 

0.97 

0.90 

0.11 

deg 

*1 

0.00 

0.00 

0.02 

deg 
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FIG. 1. AFE RAKED-CONE CONFIGURATION. 
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FIG. 4B . RFE TRAJECTORIES . 
LATITUDE. 
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FIG. 4M. AFE TRAJECTORIES , INERTIAL 

LONGITUDE OF ASCENDING NODE. 
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FIG. 4N ■ AFE TRAJECTORIES . 

INERTIAL WEDGE ANGLE. 





